“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1998-06-01 


The role of salinity in equatorial mixed layers 


Stougard, Pegeen O'Neil 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/8930 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


f (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

| | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 
























tip avonhse 






‘ ype FY ery Ld a MES pgs da 5 tat res bri > a H Te ae eal ye Ry fy gp Od tory ey 
arty tr ‘ i A dityabas sd shsbede are . OT Re 6 Nor oR ee 
RPSL ea i ! at rrr Ae Fy F <i n PE rye a F lias iterate ‘a i$ RRS ; Ae eats 

































“* ae | Paar a : ry t] ry rey ae 
LE LL TRE Ayeeyin JE ! ae ref i 0 rT Lanaetinn tet a) ‘ t Lit a te) pr atta Epa 
Sea fia seg eal va Hate | Fy ba Rea +L on Rr BH eri Th My AE ot | y ae Jad, ie rey re] ree tener Sage Seer ies 
ST oa ae ra il } Ped RTs Rid 4 ae a I. yi a 4 retpe) Reb born Sted 
ee et pein dre Arata fy) Ray ar 5 arr ert ety 
REE ey RPM Neste; Par net ies ie re eat TO ‘ terbekeeal i eg Ae rte rat ‘ q soe wees aes nt 
A i . ait if Paes Poa ray ‘ a c Pal el ’ ot Cho 
| ed ey Ry cara aaa Tea le nasty 


me ee! Lee it Fttthee 


ve rs aosased Oy rl 


Ory CU a 


Sake Py tC re BM a ok erste PAR ull, 
a, : a v mer rf i eas ty be pate “Zl beth 


r LY 
He H Thee | bs ha mH t) ry A a P OF : re an ete 

Se ~6=©6MNPS ARCHIVE Se eee re ei ee ee erie nee tig. 

ae Pico ae Co oe ‘ H } F 4 i hl oy ; fe er oR wrt Ra 
fy oe “ ; f nt Bs na ibs Tah ; a ehid z Tot ee 

i 1998.06 Tra near nee on ence 

hier) o : ie ara r reer] 
rps 


“ * b F Ne & H Le p eet pe er rae 
cea STOUGARD, P. Sauls a rahateasta ae ee aH es io sue ie 








Uy ey pra ar » ty bed n 
na re ‘ Poe ayer St eu AHR Raa A > re Hr 
aM LOR rar tras 3 g Pat 


ie 7 ry 5 i ard E 5 P wt “er. 
Saat ad f) pana rpbasiarerd ‘4 Meee i Ther Hi a5 baxter if oa ? L ni 7m Rarer 
ony +h im dol. ett . rary F He oi Prorat a . z wa Yat) Pa pes st Sea ren maT 













Nice! 
Reon e Bi Pah 
































F Cateye s reels eH Pn) O ne eer) 
(ote ey Te SE re eat ctr a 14 Lava bes, wine Leta teed aoe ah fg nvet Hat ice Ty P f 7 * phe 
ee rer CASE RO Vb te at hail ies rhe Pov U wd a 44 HE no Pedlad evita teout PtP ’ t 4 ir] ay Tht ed a1. cy rare i pyrerh va 
baa) retiree ae Be ey 5 a RCT Pa ie ne PPT Eee tes Wee iy ieee ‘ Hy Mths ba wtalog ited +) oy PPT ee ey ee 
eat eee Tee oe | t ri ue oy 4 L 5 
erictinves Rettig aCe Nears yt hae wbate i Le yy iy Pie 4 SL ed ot . i Pt 
rte tk Maan Pe ed 4 u Cd 4 





shies | 
ish 





on | rrr SPT) ae oe ae ee 
eek CL ere eed Pre eC) ey Ei a 
1 VqNhow bad Lake terre 
art Beth pape se 5 
TRiege Tavs eeesd 


2 
ee ee he ort 4 
a ra re 
eC tation pete 










rs 
. Hi 
oy) a ey. i wr) 
are Vea Pt heuyeobere 
gittw ose beage Uefa 
.) Py) Ce Pare * 











































i “Ana rr) 
A a a ar 

Ste’ “ee y 

a i | 0 i] 

7 a. A Pi Pord art it ae 

SHIH aT Pree wena he oh TE rf UL Be en ary rales ees re erect srs hart viet a GE m ; re ae ae ype — ote 

Tare a ae Br eyed : Ce Ur riet is Yr per Mt by ate jas nade 4) H int ey F re a ca ra st 

ie a ae Clr Uh id 4 F oT TE F eo as pot ae ipa Vaiekae vin PAM et Ye Hoe me? Yi Crean ; Pasa eet AP re oe apace poh 


re rt] wa 


44 rns Pi ry Sere 
ony Ike es) eretrirey st Lp eed ee prem Lead teat se ~ 
r 





























ih ay ee 
Viren aecry Ne Me He siytg tober aah 98 atte By 5 he ot od ey ene - 

rary } rn wal Tip rays fi it + nN PAST er or | ] reve: ry be 

Ter ie eA Teen Teeny rey its er , He erry | et sy 
ret Mere he tt TR Meee Vato oe Hore aA TEAL Peta by : " r i . Con BEA en oes ids ‘ abe ts i a 5 ; ve. r 

VEERVCM WOT, PRC La SPPrar Pir yea aes Pate a Sete PEL Ser ae PRPC PL UTERYT ESET Wet RL PENS ay aot aIubargtybacgtaderat Pan Hoe oP “ : a] bi aera + ae 5 te 
PUT TPT th MRR OMv a ELE Set tay eee tele ett my oe oe toe, ret vet ry ore] eT Hi Meee arin Md “6 . “i ry ed a a ye . 
OS tia ae oe ek es ahh pen Ror ee ee 0 Oe Pere WS on he MerUrer eure) He a -¥h ‘bres s tied pcs 
eT Mane Ca IGN VOC ees CELI BRT OTT Be re Brtrrcah oh i tray i Me Vee ir ry LPL) 
ad ee ee ey fame sig a Fel’ PAC ORP ae Bi Mt are i be Cr ear hs 


St iie Eon UL Cer ea) Ate sr RUAN Feu a a ee ve PTR EAS, 
Us S ptl uae rele reas eee a PA ae ay 
Ore eT evar Tie Pied | Ae Arh pivte OR 0 
SS a ae et ite tear Le a Nir ea S| 

7 eS ener rey 4 AU eee hs 

TY Se Pate Ca zee ey CY Ok 

oP ooh Bd Dee Oats Tara 








rere <) a 
or ea eee oY cr 


























Phe he RY Lo Lee 


































a phakl a ee teeey oa eh ae | ie PE oe a 
Ee oe pleas HAT CPP Ul Peat hes Sine ey Bo tig! A tt { : . i f 3 . 
NAPE TENSE Taree TELL rary aah ret aS bi ; har ert Sraek ‘ or Z KS o8 
rime it Ey Tse qa 4 Pail a} aT Te Le oe LR Her , Sata 
reo eee rat Pee ar by " rahe ae ig Me, " A 
nt ‘ 5 4 Py f ¢h A im . 
yl re Tease onda 
RHA eck error Tena 3 con 
rh ape adds id aage Oh Haan nee aero a be “4 | : teat . y Ta oocee 
a ‘ r aed Wedges PPegreric t) et by = ro wea 
a yt de sar ue “ ay al MRE ae | ry ry] 71 Ree Gy Poet Pern 4 | Yaea . ry 5 rr Lefivioledo tee: 
Sbea TP ae ea DP PS Sap de tg ya ; , e 4 ae 
ce igh ee? eee be 





sf a 
= 7 nee A eh hae ler petigh res 
Howiretatgelerse ts’ Of fags. 
E rr " 4 vad WRada® bond 
B ae PET ey Cn Re Ee Me aT Sy 


> a 
eee ae Ashe 
is rhs 











ne 







































i 
vy Sos Ory oe ee Pe A) + $44 marries 

A Sal! tee ta | ara YH carer 

tre ees eee Malou. ties 3] a) 
i ee era PLY Pea kas 

gdoQedrta lite g ” 4 thas a 

Lt. Po ay ery ‘Tre Ore PLE er oar fF ry LE ate 

Sy CT ere ie Oe sd be, ee Te | oe aryagShs ren 

nefacrdts teen’ rr) Peers Pe ‘ ke | "aS AP start 

fe A wt) Moet epee | ‘ 











My 3 
ao 
Pat nee Pr uy ra 


Ee wD 






4 













: é pit 
7 ett he] Wiel whgic Parte eo 
. Ee 




































ire Come eh Bh eye 2 d cSt A(ade 5 oe Tad 

SY aS Pe ep i ie 4) Arb ae @ PPR Toren Pepa tye Pan eee . 
A TARRES eee En ch area Pures i) tu re RR ah uly an eT ey ‘ 

a) en Pris ee ee a | ote Tdiky ee St ey a a ee Wi et H o% aS at, 

t e where “Ugts fvin picdse! *G.i or) bed Pl art realy ie Ah oor na) ‘ 
Pan aCe ( $a gh 44g 08 pores 4 ne ar od peg de eres 
Lae Te ae ad 2 Pa le eb ot ub AOL Shay be oe Bo ee oe 

fr ALe tr at 4 A Th! “ANE ay ra hy “fsena. wie varaty Pee Ay tae a 

rt some ke Step es re et ad i] 
H 


rd 


eH 
ah ake Lae 


Pook bed 


Het Male * A 








































































re ea 
J » 
ar 
7 £ 
ta 
| H 
Po , 
rn 
arn Ae ari cr 
MTL 7 arr irk 
’ 142 fF has 
He 
LS 
aft 3 . 
r rt 
4 ‘ 
7 
at, sl 
oF i 
! ‘3 
Ca ‘ $e 
‘ ce 
Ae rs 
ws t 
A i: 
‘ , hs 
; aT 
P 
a tegee ; 
ore 4 
Le 3 
SG tt pd 
att B bi 
4 : Y ¢ - 1 fi 
J 1 ey . cs ee iS Tr a) ace lee sy ws 1 hy 
. ; q a fe eT ‘al oe) rE 
i a Lat i r Pa hk aren ib F AA 4: a 4" yh AK i ey 
Ch Pa Cs Le Bi f: 
Arkh ete ete ah abating g deem F top Ch ; rat ts Bey : 21 
hea Wee Sey SRLS Ty aes OL 0Y | erate : 4 qd4yy ie ea es a r : : ” he \ CMON ITE aaa t mri wos 
ae Loh OL ee Vy Sanh a hd ae : He $ TF ; 4 Fs i 7 : ir teh i = Jom ieg-asqutart Am Daeg 
et ri oi Pperien: Pek Wed Om ae Pai? rE iy P ; Y ; ; atte’ : 
Py t- Pa A Pe et er A ee 
a5 MN %.¢ B S y ie " i ‘ , . if ¥ 
ff " c Ant , Hr Me | 4 Fhe " a EL in bt eae ohag! ; f “fi 2 rit prot ¢ Pe oe , i pve 
, H " A ren ie) er re r A 4 en 5 UH : Pat ipa tra) 4 , J ee, By ert ei i) ch! : ye 7 oe ahi 
Uh Te ‘ i : ’ ; $; A ee : 5 a H ELC Gens are fe “t i is 
iP VOU Is if 5 4 i 4 ce ‘ a ue ‘ 3 4 oe J 5 ree n 1 b aes oy 
a ‘ aD A 4 oe 7 hi ‘ « 5 * R ’ k : bat i . nA f PL be ee 
veoee P i ; ; : : Sear ey ey i 
m 5" i au 7 a i é , whe a, f 7 2 Te ' } 
Daca Ure O ES Oe Ia 2 er ae f : che 3 k Tye Ge aro aa f Leer 
veur ey a 7 A eee TT ¥ b } 
Ca eres? ae 
| ' a 
TT vent 
i r te? na yi Tare ths 
Fan Y Ln WL eee a Le ee 0 ay! a ee iG 
‘pp a) dl B a et ete 1% i iter ‘ oe Be Oe rh DER) ; c A F re ea . ( ' . he , . Aa : “4 
are ae Rare. 4) a; ae oe | arty rhe tt} Ut Aas Pays) Ps rey ts } ' Tae te! . 7 oe > at) J dade Weal bok edhe 
fadnge fit J’ Fore es ae ae ba eT) eh) eh ad] ah J be Lae yas ay : ma AIP a es a . oe ea P 7 " t ae 1h Lite ald tt oh hte 
avo yl PP RR TRY fo Lie bk RP uk tet br a8 4 vey Pepe pir, ee 4 ‘AS ? ay 4 243 Paso & bee bend bat 4 cas tlh fe ia nti 
‘ ee pois ee ig marie i may AL Pee SO aE eed riers ; p i Agr goo tie: et 5 Reb bse! eee 
a! 4! Brie ae Th : 1 t ¢ Ap Le iH ous sPa\ se? Lees hao 
git an et a et th BTU HO Pent aL ee Pct Pe BH oe ARE MN Phe Pesci ied Sal. 
; cs a oN 5 Lr Pit he: : : i, : a - Ry j % o é i A u ree + Pont i a % K Kein Pe boas 4 bh off hes 2 Meh f 
t rae 3s ‘ PF i ‘ Bh f C Bret PA ak 
, Dv a ie Pett oi Fake 2 
RAL Ai eh BLD aL Bote a fost ighcae hate aE 
Lia ae Lee Db gered gs coast * 4 


EN cg eat raed hJ 
eh dele ett aed ple dad 
VENAD PTE OT 
iO aa? i 38 
Fee q oe Tiare Rhee erke et bys 
J iy iF , Pa nae F Fo Yj b Me - To > oa baa 
ae ie 
aus i] t i) fa ‘oe? yee Py 
h seH9 rh ha] Daal F Ci ee ied bate ER Lt ha td be lk bie hats ot 






fa . 1d Les tee ah 
Ta EN 1 ith ate bela hs Ey airs | Br = a 4 ; MA peed nt Eat ‘ezcrmeny 
fee yar Pha 2 2) ty DT Ptee Be patecta De hs bh A atl ad 


Nae a ie ae E t 2 " kd be ea g a 

PLT i S F Ate teh Date Wh at hee 

ere SET RRO RE MERC eat abrairt a pater et 
7 re Psy ed od eh co + t ; ein tet Cte dee 

feeree vite tay hy Bu 4* a Setar patetes: Freese 

Me tsa yearg : 













PT 


BY 
i Ate yd Eee ay 


Pe 
op en 


ry 
er 
rd 
~ 
we 
a 
eed 









» A 4 f at + ee We Daal s fens 

<4 et ” L it > aie) 5 Ba ‘ rea be foamed 

Sa be H he a x , , Hy i ‘ A f +s PSE i Ea 
rt ha beg EP : Les ae ry 


Ws Ha ee 
3 se F 
oD a tee 










Heelan 
























es tn rie AT P ieee d xt od 
aa a re err ytery jor daes 
rae) ph Gee edd dba pede Sopp eet a et Sr aes be Pir Mie : fe > he iS - 
farb vi ecs f mata a eh ane ed Yad resus ‘ Ale ; r bis . ” ; BP im bs hk Dae ved 
= - es « OER thee f ; mag] rTiee & = bs ; O 
; a eseelaen cael HF i? a SORRGATT Be #4 Ta Pre t r ry bec ENP k eds EL Daa bar cobs he 96- oee iaa nese 






42) ste bh Mthy Hey Ley a 





tr eel Re . PPC HbaeE 
































Ded hs ad Fh Pie Wt bd, dy copia adh belehanth tp died 
Ase : r vee rn) Pee Ly Los a ha Deak) sana Me 
art ae | py LAE a H A ? ae) ot . Wins PH Fa ae 4 =H Rad ad TwTe Acad a ee Ct is oe 
: inden eeaueat Rett oe Peete For erg 4ty Ss j ee SRE eh he! tbh Refined 
n fi ar ta ‘4 Boa Bh if o 
eave Perrier sce Chath: | 
5. , ET oe : 


fhe i 

ert ert ts 
eRe eats he beet = 
eee LL dk inlaid 





- re] BS Rane rr or) a f' 4 vintar miricsnenk es 
witeed APeeY ¢ ‘ i rE PAaEd 2 ot, id Arti nha 
é by rd z hud! 
4 FH oa) 
a te he Pie® ,F 5 | ip, | pe ae 
i: ‘Le rH he * Deepen ass Bara E are 
Pris} phe tay esate Ceo, 
HA i 4 1 A rer ae CL 4 ite eS 
re ars Bt AP Pe ete) a te ri he ie ie St Phe 4 LT Lop pA 
pr UTS a ta Tae age. a SD , oy he H ia Ue Bt orn: f Ml Py ry nt B He : f ere 
1 et Pee See oe ike ORE : De LR AD Ne eee Ha A pte Be ae Re cites Sichiecaier tts i i rio ey * cee 
bu tatoe J a MU wet AOL 7 ‘ ‘ J wha 9 bs Fs oR rer at : - 
] v Da LS Bea! by ho oty Mins , ry be s ¢ eee nS val pao At 5 sy tee ee SS len 
ar hh Wee PEC reed dP ay ri ae Lae y A , ier t e 4 4 p 4 4 te a Sh v oy ae = 
Pr ee Fl Tees i wh Retreh Spas Sete ok oy 


f 4 , r ere r ' : P t F k > 
0 qi 7 4 r ; 4 ir l ; i} H 
pPyotsh, 6 fd 8 by Mitweegintatet yt ws » " re Pe! ts 3) de Fig 1 Cat, NO dey Sa rngt gs PEs aera ht sir] f at ‘ PF i tna ee blk Shek 


DUDLEY KNOX LIBRARY 
ME&VAL POSTGRADUATE SCHOOL 
UITEREY 2S 33943-5101 











NAVAL POSTGRADUATE SCHOOL 
Monterey, California 





THESIS 


THE ROLE OF SALINITY IN EQUATORIAL 
MIXED LAYERS 


by 
Pegeen O’Neil Stougard 


June 1998 


| Thesis Advisor: Roland W. Garwood, Jr. 
Co-Advisor: Arlene A. Guest 





Approved for public release; distribution is unlimited. 





—= = ae — SS NN a — ——— _— 


| REPORT DOCUMENTATION PAGE jem aemoesonne omtase | 


H Public reporting burden for this collecton of informaton 1s esumated to average | hour per response, including the tinie for reviewing instruction, searching exisung data 
sources, gathering and mamtaining the data needed, and completing and reviewing the collecton of information. Send comments regarding this burden estimate or any other 


aspect of ts collection of micrmativiz, mchidnig suggestions tor reduciig this burder, t6 Wastmigton readyua ters Services, birectuvaic for Ldurmatun Operatwus and f 
Reports, 1215 Jefferson Dams Highway, Surte 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reducton Project (0704-0188) 
Washmeton DC 20503. 


1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 
June 1998 Master’s Thesis | 
TITLE AND SUBTITLE THE ROLE OF SALINITY IN EQUATORIAL FUNDING NUMBERS 
aaa LAYERS, UNCLASSIFIED 


AUTHOR(S) 6. AUTHOR(S)PegeenO’NeilStousard O’Neil Stougard 


PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) . PERFORMING 
Naval Postgraduate School ORGANIZATION 
Monterey CA 93943-5000 REPORT NUMBER 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


11. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the 
official policy or position of the Department of Defense or the U.S. Government. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 
Approved for public release; distribution is unlimited. 


13. ABSTRACT (maximum 200 words) 

The purpose of this study was to understand the role of surface salinity flux in changing heat exchange 
between the ocean and the atmosphere by means of its effect on mixed layer dynamics. This was 
accomplished by a series of thirty-day mixed layer experiments using the one-dimensional Naval 
Postgraduate School (NPS) mixed layer model. Results from the NPS mixed layer model, forced with 
both idealized and in situ data from the western equatorial Pacific Ocean, demonstrated that salinity can 
play a significant role in potentially changing the surface heat flux, with its effect on the mixed layer 
depth and mixed layer temperature. Precipitation stabilized the mixed layer by creating a barrier layer, 












which slowed entrainment. The net accumulation of rain was found to be an important source of 
buoyancy that reduces entrainment by subsequent wind mixing events. 





14. SUBJECT TERMS 15. NUMBER OF 
Oceanic Mixed Layer, Salinity, Ocean Models PAGES 67 
17. SECURITY CLASSIFICA- | 18. SECURITY CLASSIFI- 19. SECURITY CLASSIFICA- | 20. LIMITATION OF 
TION OF REPORT CATION OF THIS PAGE TION OF ABSTRACT ABSTRACT 
Unclassified Unclassified Unclassified UE 
NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 


Prescribed by ANSI Std, 239-18 298-102 


11 


Approved for public release; distribution is unlimited. 
THE ROLE OF SALINITY IN EQUATORIAL MIXED LAYERS 


Pegeen O’Neil Stougard 
; 
Lieutenant Commander, United States Navy 
B.S., United States Naval Academy, 1988 


Submitted in partial fulfillment 
of the requirements for the degree of 


MASTER OF SCIENCE IN METEOROLOGY AND PHYSICAL 
OCEANOGRAPHY 


from the 


NAVAL POSTGRADUATE SCHOOL 
June 1998 





DUDLEY KNOX LIBRARY 
NAVA). PCS TGRADUATE SCHOOL 
ABSTRACT MONTEREY CA 93943-5104 
The purpose of this study was to understand the role of surface salinity flux in changing 

heat exchange between the ocean and the atmosphere by means of its effect on mixed layer | 
dynamics. This was accomplished by a series of thirty-day mixed layer experiments using the 
one-dimensional Naval Postgraduate School (NPS) mixed layer model. Results from the NPS 
mixed layer model, forced with both idealized and in situ data from the western equatorial 
Pacific Ocean, demonstrated that salinity can play a significant role in potentially changing the 
surface heat flux, with its effect on the mixed layer depth and mixed layer temperature. 
Precipitation stabilized the mixed layer by creating a barrier layer, which slowed entrainment. 


The net accumulation of rain was found to be an important source of buoyancy that reduces 


entrainment by subsequent wind mixing events. 


TABLE OF CONTENTS 


1 NER OI UC PIO IN eae i eee eee oe eee ] 
A TO Ae Ca oy te era ees ete) gee are se a RR gee l 
B. RITERA TORE RE VIEW owe ee ee ee ee aa 3 

HINES 1 -DINIODEy es ee ee ek ee 7 
A. MO DED) Visi TE We ecg ieee eee eee eee i 
B. DA VA SOUR@ES coe yee eee eee ee ee 10 

I. SALINITY EGEECTS ON NES*M@QDEL go 242 ee es aw 13 
A. BASIC SATIN Vee BEBCES oy 6456 se ee ee eet ss eee 15 

1. Constant Salinity (Case 1) vs. Climatology Salinity Profile (Case 
OR lean es eee eee Ge a, 15 

2 NWyind. Biirects (Gases 3 anar4) 22 2 es ee ee ee ee 20 
B. PRECIPELA TION PRREGTS «0.654542 64e- 64 649 +e eae 
1. Steady Precipitation (ase 5)... 2. ss 34445. see 23 
2. Diumnial Radiative Forcing . 9... 2.22.2... 1.2 tee 23 
a. Four Hour Period of Precipitation Starting at Local Noon 
(Maximum Heating) (Case 6) ................. ZS 

b. Four Hour Period of Precipitation Starting at Local 
Midnight (No Heating) (Case 7) ............... 25 

C. Four Hour Period of Precipitation Starting at 0600 Local 
(iMormine aim) i(Case 8) . ... seein ee 28 

d. Four Hour Period of Precipitation Starting at 1800 Local 
(Eveniic Rain) (Case ooo. oa. ee ee 28 

3. Heavy vs. Light Precipitation (Cases 6 vs 10)........... 28 
a. Fifteen Hour Period of Precipitation at Local Noon (Case 

BO ak ee a eee a See eee 28 

4. Wind Bitects 9 2) e eee 4 eee 32 
a. Light Precipitation (Case 10 vs Case 11) ......... 32 

le Heavy Precipitation (Case 6 vs Case 12) ......... 34 

c COARE FORGING: oi. 4526 2 ae 1 o eee eee 34 
l. Hull Orcine CCACE Ml 3a care ee eg ee 36 
Z Partial Forcing - No Precipitation (Case 14) ............ 38 


Vil 


3: AVCTAC CURE ORCHID ents a 0 yernt eatin. swe ie te a 
a. Averaged Wind and Precipitation Forcing Only (Case 15) 
ne ee ee 43 
b. Averaged Radiation, Wind, and Precipitation Forcing (Case 
OY oan em cua ee eee ee 46 
IV. A BRIEF LOOK AT LARGE-EDDY SIMULATION, DISCUSSION AND 
RECOMMENDATIONS wr. oa cs ote uc eee eee err 49 
List OF RBRERENGES. ....: 205) See <8 ee 53 
DNS ETAT STS ERS CL) INOS a geass ese racer eee cane rg 55 


ACKNOWLEDGMENT 


The computational costs of this thesis were supported by NOAA and NSF Physical 


Oceanography under Award Number 9413292. 


I. INTRODUCTION 


The purpose of this study is to understand the role of surface salinity flux in 
changing heat exchange between the ocean and the atmosphere by means of its effect on 
mixed layer dynamics. This study will be accomplished by a series of thirty-day mixed 
layer experiments using the one-dimensional Naval Postgraduate School (NPS) mixed 
layer model, which will be discussed in full later. The hypothesis being tested 1s that 
without significant changes in the surface heat flux, precipitation can by itself vertically 
redistribute the ocean’s heat content by salinity’s effect upon stability. 

The primary area of focus is the western equatorial Pacific’s warm pool, which 
was the focus of Tropical Ocean Global Atmosphere (TOGA) Coupled Ocean- 
Atmosphere Response Experiment (COARE). This study is part of the TOGA COARE 
research program that was funded by the National Science Foundation (NSF) and the 
National Oceanic and Atmospheric Administration (NOAA). 

A. TOGA COARE 

TOGA was a ten year project (1985-1994) whose major objective was to develop 
an Ocean observing system to support studies of large scale ocean-atmosphere interactions 
on seasonal to interannual time scales. This ocean observing system was completed in 


December 1994 and is known as the Tropical Atmospheric Ocean (TAQ) array (Figure 1). 


The array, which spans the equatorial Pacific, consists of nearly seventy moored buoys 
capable of measuring oceanographic and surface meteorological data. (McPhaden, 1995) 
This thesis will focus on two of the four main goals of COARE. These goals are 


to describe and understand (1) the main processes for coupling the western Pacific 
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Figure 1. TAO Array. (http://www.pmel.noaa.gov/toga-tao/gif/mooring.gif) 


Ocean’s warm pool and atmosphere and,(2) the ocean’s response to buoyancy and wind- 
stress forcing in the same area. Achieving these goals are important for understanding the 
western equatorial Pacific since it is the location of some of the warmest sea surface 
temperatures in the open ocean, and it has the largest annual precipitation and latent heat 


release (Webster and Lukas, 1992). 


B. LITERATURE REVIEW 

How does salinity affect the mixed layer? Miller (1976) was the first to include 
salinity in a mixed layer model. He used the Kraus and Turner (1967) model which had 
been modified by Denman (1973) to include dissipation. Miller added the conservation of 
salt to the equations to examine salinity’s effects on the mixed layer. His results showed 
that a homogeneous saline layer significantly affects the depth and temperature of the 
mixed layer. The increased difference between the mixed layer’s temperature and salinity 
and the underlying temperature and salinity cause an increased density gradient. This 
increased density gradient leads to a change in heating and cooling characteristics because 
it slows the deepening of the mixed layer. Including salinity also tends to reduce the 
cooling of the mixed layer because of entrainment across the mixed layer boundary. With 
a salinity difference at the lower boundary and a significant net cooling at the air-sea 
interface, it is possible to have a mixed layer that is cooler than the layer below. It will not 
be unstable and overturn as it would if only temperature were contributing to the 
buoyancy. Depending on which dominates, surface flux or entrainment flux, adding 
salinity can cause an increase or decrease in the cooling of the mixed layer. (Miller, 1976) 

Why has salinity been overlooked in equatorial Pacific ocean modeling even after 
its effects on the ocean in general were shown? There are two main reasons. Firstly, due 
to poor data coverage in the Pacific, evaporation and precipitation over the tropical ocean 
are not known with any degree of accuracy. Secondly, because the thermal expansion 
coefficient reaches its maximum value at the equator, temperature was thought to have the 


dominant effect on the dynamics. Salinity does not have a-direct effect on the atmosphere 


and was assumed to have a negligible effect on tropical dynamics (Cooper, 1988). 

Cooper (1988) showed that complete density data, comprised of both temperature 
and salinity information, caused the basic tropical ocean circulation model to produce 
accurate results during a sixty day simulation. Using just one of the density parameters, 
either temperature or salinity, gave poorer results than if no density data had been 
available. The differences due to the individual temperature and salinity gradients are 
larger than the combined density gradient because of the relationship between temperature 
and salinity in the equation of state. Cooper also showed that the modeled temperature 
and surface velocity errors increase when salinity gradients are not included. Temperature 
differs by a maximum amount of 2°C and the surface velocities differ by up to 40 cm/s 
without salinity in the circulation models. 

In the past, mixed layer and thermocline depths have been determined from data by 
using a temperature gradient criterion or a net temperature decrease from the surface 
temperature. Levitus (1982) used this temperature gradient method to define the mixed 
layer depth. Knowing the importance of density differences, he also employed a change of 
density criterion. This density criterion was used because of the stabilizing effect of 
Salinity in the upper ocean. Lukas and Lindstrom (1991) used the Levitus method in their 
study of the western equatorial Pacific Ocean mixed layer. 

Lukas and Lindstrom (1991), using CTD measurements from two cruises, found 
an average mixed layer depth of 29 m. In earlier studies, the western equatorial Pacific 
Ocean mixed layer was calculated to be about 100 m, three times deeper. The earlier 


calculations were so different because they were only using temperature gradients as a 


criterion and the effects of salinity were not included in their calculations. Using a 
histogram of the differences in the thickness of the isohaline layer and the isothermal layer, 
Lukas and Lindstrom found that more than fifty percent of the vertical profiles had 
shallower isohaline layers. This finding suggests that precipitation effects were important 
in determining the mixed layer depth. Precipitation formed a stable, less dense layer of 
surface water referred to as a “barrier layer.” This barner layer plays a major role in 
inhibiting entrainment. They also found that the shallower observed mixed layer depth 
appears to result from a strong positive buoyancy forcing associated with heavy 
precipitation combined with intermittent wind forcing in the equatorial Pacific Ocean. 

Miller (1976) also examined precipitation events and their effects on the mixed 
layer in his initial salinity study. He found that the depth of the new mixed layer caused by 
precipitation depends greatly on the amount of wind stirring during the event. Higher 
winds produce a deeper mixed layer than lower wind speeds when they happened during a 
precipitation event. The resultant mixed layer depth depends on the duration of the 
precipitation and the wind events. Finally Miller found that when heavy precipitation 
occurs over the tropical ocean, a stable mixed layer is developed at the ocean’s surface. In 
addition, Flament and Sawyer (1995) noted that the contribution of a 1 cm/day 


precipitation rate to the density structure is equivalent to a heat flux of 40 W/m’. 


II. NPS 1-D MODEL 


A. MODEL OVERVIEW 

The NPS one-dimensional mixed layer model ts a revised version of the Garwood 
(1977) model. Garwood’s entrainment model considers both turbulent erosion and 
dynamic instability. The version here used has been modified to include salinity. Equation 
(1) is the basic definition for the rate of change of the mixed layer depth (h), and its 


solution is fundamental in the model. 


oh = 

—-=w,-W_,, 

Or e s=-h (1) 
where Wi , 1s upwelling/downwelling. Assuming no mean vertical motion, W=0, the 


equation becomes: 


oh _,, 2 
Ot ©) 


where w, is entrainment velocity. The entrainment velocity is always greater than or 
equal to zero because the water column cannot be unmixed. Defining total turbulent 
kinetic energy (TKE), 
E=u2?+y2+w® (3) 

the vertically-integrated total TKE equation is: 

< (HE) = 2mu,> + (uw +v)w, - ghAbw, - ghbw'| aD) (4) 
The terms from left to right are the storage or time rate of change of net TKE, wind stress 
shear production, entrainment shear production, entrainment buoyancy damping, surface 


buoyancy damping or production, and viscous dissipation. The dimensionless wind- 


stirring constant m, is estimated to be 6. The friction velocity, u, is equal to /t/p, where 
wind stress, t, equals p.C,|u,)|° and p,, Cp. and u,, are the density of air, the drag 
coefficient, and the wind at 10 m above the ocean surface, respectively. The variables u 
and v are wind driven currents and are assumed to be equal to zero below the mixed layer. 
Gravity is represented by the constant g. Net dissipation is defined as 


D = 2m fee where m,=1. 


Changes in salinity cause density differences that in turn change the buoyancy flux. 
Buoyancy is negatively correlated with density (p), b=g(p,- P)/P, where p, is a constant 
representative density. Thus, a more buoyant particle has less weight. Temperature and 
salinity are positively correlated so that warmer water tends to be saltier. Usually, salinity 
reduces the density gradients because it counteracts the thermal effects on density. This is 
easily seen in the equation of state, p=p,[1-a(7-T,)+B(S-S,)]. In the equation p, is the 
density of sea water at a temperature T., and salinity S, . The thermal expansion 
coefficient is w, and B is the salinity expansion coefficient. The values of T, and S, are 
representative values of temperature and salinity for the case being studied. The thermal 
expansion coefficient is a function of latitude. In the midlatitudes, « is approximately 0.2 
and increases to 0.3 near the equator. Therefore, in a tropical ocean one psu salinity 
change causes a density change equivalent to that caused by a 2.5°C temperature change. 

Dividing the NPS model’s buoyancy (b) term into its mixed layer temperature (T) 
and salinity (S) constituents gives the buoyancy jump at the lower boundary of the mixed 


layer, 


Ab = (agAT-BgAS) , ag=0.3cm/s/C and Bg=0.75 cm/s/psu (5) 


where 





N= 7 andwkS = S-S_ (6) 
The surface buoyancy flux is 
b/w(0) = agT'w‘, - BgS'w’, (7) 
where 
Pg Oe fond ¢C 
T'w'y = ac: S'w'y = (Pr-Ev)S and Q,,,=Q3+O,+9O,-Q; (8) 


Q,,;18 net radiation, p is ocean density (approximated as 1.028 g/cm’), C,, is specific heat, 
Pr is precipitation rate, Ev is evaporation rate, Q, is back radiation, Q,, is sensible heat 
flux, Q, is latent heat flux, and Q, is downward solar radiation. 


Solving for w, using the quasi-steady state mode, where the total TKE changes 


with time but the storage terms are negligible, 


ae 

—(hE) = 0. 

a E) (9) 
Replacing buoyancy with temperature and salinity by substituting Equations 5 through 8 


into Equation 4, the result is: 


9, 


2mu,> - gh{a—-B(Pr-Ev)S] - D Az. when 2@>0 
Aw, = See a (10) 
gh(aAT-BAS)-(u" +v") A=0 when a 
t 


The primary difference between Equation 10 and the equation Miller (1976) used is the 


explicit inclusion of dissipation, D. The modified Denman (1973) equations, which 


originated from Kraus and Turner (1967), accounted for dissipation by incorporating its 
effects into the constant used in the wind stress production portion of the equation. This 
incorporation reduced the wind effects, Le., 2m,u> -D=mu? in Miller. This 
approximation is good only if D is always directly proportional to u>. IfD depends upon 
buoyancy flux, then Miller’s parameterization may have significant error. 

The NPS mixed layer model is initialized with a temperature and salinity profile 
chosen by the user. The model can have either a constant surface forcing, or the forcing 
may change at each time step of one hour via a data file. 

B. DATA SOURCES 

TOGA COARE’s Intensive Observing Period (IOP) was the primary source of 
data for this study. The area of focus for the TOGA COARE IOP was the western warm 
pool of the Pacific Ocean, located in the region contained by 10°N, 10°S, 140°E, and the 
international dateline. The IOP took place from November 1992 to February 1993. The 
center of the [OP’s data collection efforts was the Intensive Flux Array (IFA) located at 
2°S and 156°E. (Webster and Lukas, 1992) 

A temperature profile from 1.75°S and 156°E on 5 December 1992 was obtained 
from the TAO Array data and was used to initialize and compare results from the NPS 
one-dimensional model. Since there was no corresponding salinity profile to accompany 
the temperature profile from this source, another data source had to be used for the 
salinity profile. An annual average salinity profile located at 1.5°S and 156°E from 
Levitus (1982) Climatology data was used to initialize the model. The temperature and 


salinity profiles were linearly interpolated from the surface to 200 m depth to a standard 


10 


spacing of one meter. Linear interpolation was used to replace any bad data points found 
during the period of study. 

Forcing data were acquired from the Woods Hole Oceanographic Institution 
(WHOI) surface (IMET) mooring, which was part of the IFA deployed during the TOGA 
COARE JOP. The IMET buoy was located at 1.75°S and 156°E and measured net 
shortwave radiation, latent heat flux, sensible heat flux, net longwave radiation, wind 
stress, rainfall rate, and the sea surface temperature. The measured data from the IMET 


buoy was averaged to an hourly format. 


I 


Il. SALINITY EFFECTS ON NPS MODEL 


Using the NPS 1-D mixed layer model, various forcing combinations, listed in 
Table 1, will be tested. The first two sections of this chapter comprise a sensitivity study 
of basic salinity effects and the effects of precipitation. These forcing combinations are 
not normally isolated in nature. The various forcing parameters have been chosen so their 
individual and combined effects on the mixed layer can be systematically examined. The 
first section isolates the effects of a salinity gradient with negligible winds in cases | and 2 
and with moderately strong winds in cases 3 and 4. With just salinity gradients affecting 
Equation 10, a reduction of mixed layer temperatures should be evident. The second 
section focuses on precipitation effects under various conditions, including timing of the 
precipitation, changes in the precipitation rate, and light vs moderately strong winds. With 
the inclusion of excess precipitation, Pr-Ev>0, the entrainment rate should be reduced due 
to wind stress and Q, alone. Precipitation will have a short term effect on the mixed layer 
through the surface buoyancy term in Equation 10. It also has a long term effect on the 
emits rate via the entrainment buoyancy term. This long term effect is due to the 
precipitation induced change in the salinity gradient, which results in increased 
Stratification in the upper ocean. Wind mixing will have an increased effect when salinity 
is added because of the reduction of the surface buoyancy flux term. But with the 
inclusion of salinity, by either a positive salinity gradient or changes caused by excess 
precipitation, the wind mixing effect will be reduced because of the influence of the 


increased salinity gradient at the entrainment zone. The last section will show the different 
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Table 1. Forcing used for NPS one-dimensional model. 
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aspects of in situ forcing and their effects on the mixed layer. This section includes 
comparing precipitation vs no precipitation and instantaneous vs averaged forcing. 
A. BASIC SALINITY EFFECTS 

The model is initialized with a temperature profile from 1.75°S and 156°E on 5 
December via a bathythermograph (BT) and either a constant salinity profile to simulate 
no salinity in the model or a Levitus climatology salinity profile from 1.5°S and 156°E 
(Figure 2). 

i Constant Salinity (Case 1) vs. Climatology Salinity Profile (Case 2) 

Using no precipitation and a minimal wind of t = 0.01 dynes/cm?, the 


precipitation/evaporation portion of Equation 10 becomes zero and the equation becomes: 


QO 
2m,u,-gh(a—) - D 
OS (11) 
h(agAT-BgAS) 


As u, becomes very small, the numerator of Equation 11 becomes negative. This would 
indicate a negative entrainment velocity, which is not possible. Therefore, one must 
assume no entrainment velocity by definition, see Equation 10. Setting A=0, the left-hand 
side of Equation 11 equals zero. Solving for h, the governing equation becomes: 


2m,u,” = 


OQ» (12) 


This diagnostic equation for the mixed layer depth, h, is governed by dissipation and 


surface buoyancy damping/production when wind and precipitation and evaporation are 
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Figure 2. A) December 5, 1992 BT temperature profile from 1.75°S and 156°E. B) 
Climatology Salinity Profile from 1.5°S and 156°E. C) Constant 35 psu salinity profile. 


assumed negligible. 

The mixed layer temperatures are identical in both the constant (Figure 3) and 
climatology salinity (Figure 4) cases. The mixed layer temperature oscillates with the 
diurnal heating/cooling cycle, with an increased temperature occurring right after the 
noontime maximum heating rate and decreasing temperature trend during the duration of 
the 30 day period. For the climatology case, the mixed layer salinity increases throughout 
the 30 day period. This increase is caused by the entrainment of the higher salinity water 
into the mixed layer from below. As the description implies, the constant salinity profile 
remains constant, and only evaporation or precipitation could change S. Both cases have 
a very strong diurnal pattern in their mixed layer depths, with the layer shallowing at 
maximum heating tumes. The mixed layer deepens faster with the constant salinity profile. 
The change in the deepening rate can be accounted for by looking at the changes in 
salinity over time and Equation 11. With constant salinity, the salinity jump (AS) is always 
zero. With the climatology, AS is positive throughout. Therefore, the constant salinity 
case does not affect the rate of deepening and the positive salinity jump in the climatology 
case decreases the rate of deepening because the mixing has to work against the 


climatological stratification of the water column. 
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Figure 3. Case 1 - NPS model mun using constant salinity and nominal wind stress. (t is 
wind stress in dynes/cm’, E-P is evaporation minus precipitation in mnVhr, Q,,, & Q, are 
net radiation and incoming solar radiation, respectively, in W/m?, MLD is the mixed layer 
depth in meters, ML Temp is the mixed layer kemperaturen in degrees Celsius, and ML Sal 


is the mixed layer salinity in psu). 
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Figure 4. Case 2 - NPS model run using climatology salinity and nominal wind stress. 
(Same labeling convention as Figure 3). 


2; Wind Effects (Cases 3 and 4) 
By including the wind effect and still neglecting precipitation/evaporation effects, 


Equation 10 becomes: 





Oy 
2m,u - gh(a ae 

pc, (13) 
h(agAT-BgAS) -(u" +v") 


Ww a 


With an increase in wind, the wind stress shear production should increase and the 
entrainment shear production should also increase - both causing an increase in the 
entrainment velocity. 

The overall mixed layer temperatures for cases 3, constant salinity, (Figure 5) and 
4, climatology salinity, (Figure 6) decrease at a slightly faster rate over the 30 day period 
than in cases 1 and 2 because of the additional energy provided by the wind for mixing. 
The mixed layer temperatures for the case 3 and case 4 salinity cases are similar, but in 
case 4, the overall temperature cools at a slightly slower rate. The diurnal pattern is 
damped by the added wind so there is a smaller increase in temperature after the maximum 
diurnal heating. The climatology salinity increases at a faster rate over the 30 day period 
with the addition of winds. This is because with greater mixing high salinity waters are 
mixed into the mixed layer from below. In both cases, the diurnal cycle effect on the 
mixed layer depth is damped and overall deepening rate increases. This can also be seen 
from Equation 12. With wind stress added into the equation the heating cycle is not the 
only parameter controlling the numerator. The equation also helps explain the difference 


between cases 3 and 4's overall mixed layer temperature trend. With the increased salinity 
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Figure 5. Case 3 - NPS model run using constant salinity and 1 dyne/cm? wind stress. 
(Same labeling convention as Figure 3). 
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Figure 6. Case 4 - NPS model run using climatology salinity and 1 dyne/cm? wind stress. 
(Same labeling convention as Figure 3). 


jump in case 4 the temperature jump does not have to be as large to get the same 
deepening rate. 
B. PRECIPITATION EFFECTS 

The addition of precipitation causes jess dense fresher water in the mixed layer 
over a more dense saltier ocean which has a stabilizing effect. In all precipitation 
simulations, the climatology salinity profile is used and the precipitation rate 1s maintained 
at 4.5 cm/day. 

1. Steady Precipitation (Case 5) 

When a steady precipitation of 4.5 cm/day is added, large differences occur in the 
mixed layer depth and salinity (Figure 7) when compared to case 4 with no precipitation. 
The mixed layer temperature is unchanged with the addition of precipitation. The addition 
of precipitation adds fresh water to the ocean surface which decreases the mixed layer 
salinity and has a stabilizing effect. This stabilizing effect greatly slows the overall 
deepening of the mixed layer while increasing the diurnal difference between the daily 
shallowing and deepening. In Equation 10, the salinity jump becomes negative as fresh 
water is added to the mixed layer and the salinity flux term increases causing the buoyancy 
jump to increase and the buoyancy flux term to decrease leading to a reduction of the 
mixed layer deepening rate. 

pa Diurnal Radiative Forcing 

In this section, the phase of the precipitation input relative to the diurnal heating 


cycle is varied and the length of the precipitation event is set equal to four hours each day. 
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Figure 7. Case 5 - NPS model run using constant precipitation of 4.5 cm/day | dyne/cm? 
wind stress. (Same labeling convention as Figure 3). 


a. Four Hour Period of Precipitation Starting at Local Noon 
(Maximum Heating) (Case 6) 


In this experiment (case 6) the same amount of fresh water is added as in 
the previous experiment (case 5) but it is added over a four hour period beginning at local 
noon, the time of maximum heating. The mixed layer temperature has a very damped 
diurnal cycle and an overall cooling of approximately 1°C over the 30 day period (Figure 
8). During the four hour rain event the surface salinity has a sharp decrease. During the 
rest of the heating cycle the salinity increases smoothly due to entrainment. The mixed 
layer becomes fresher by approximately 3 psu over the 30 day period. At the beginning of 
the rain event, the mixed layer shoals to 4 m and remains there until the rain event is 
complete. When the rain event is complete, the mixed layer deepens at a constant rate 
until the next rain event. The maximum depth at the start of the period is 13.5 m and 
increases very slowly over the 30 day period to 14.5 m. 


b. Four Hour Period of Precipitation Starting at Local Midnight 
(No Heating) (Case 7) 


The mixed layer temperature (Figure 9) has a similar diurnal pattern as the 
noon rain case. The mixed layer salinity follows a similar pattern as the noon rain with the 
twenty-four hour cycle, but it starts at midnight instead of noon, in phase with the rain 
event. The overall mixed layer salinity decreases the same amount over the 30 days. The 
mixed layer shoals to 6 m and holds there during the rain event following closely to the 
local noon pattern. The exception to the noon pattern is that the maximum depth of the 


mixed layer remains near 15.2 m over the 30 day period. 
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Figure 8. Case 6 - NPS model run using daily 4 hour precipitation of 4.5 cm/day starting 
at 1200L and | dyne/cm? wind stress. (Same labeling convention as Figure 3). 
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Figure 9. Case 7 - NPS model run using daily 4 hour precipitation of 4.5 cm/day starting 
at OOOOL and | dyne/cm? wind stress. (Same labeling convention as Figure 3). 


C. Four Hour Period of Precipitation Starting at 0600 Local 
(Morning Rain) (Case 8) 


In Figure 10, the mixed layer temperature still follows the same diurnal 
cycle as the previous cases. The mixed layer salinity follows the same pattern with a shift 
in phase. The mixed layer shoals to 6 m at the start of the precipitation. When the 
daytime heating starts, it shoals further to 5 m for the rest of the rain event. The 
maximum mixed layer depth reached before the next rain event was 16 m. 


d. Four Hour Period of Precipitation Starting at 1800 Local 
(Evening Rain) (Case 9) 


In Figure 11, the mixed layer temperature is similar to the other cases. The 
mixed layer salinity has the same pattern as the other four hour rain events with a phase 
shift so the salinity decrease occurs during the 1800L rain event. The mixed layer depth is 
the same as the midnight rain event. 

By only varying the phase of the precipitation relative to the diurnal heating, no 
major changes in the mixed layer were seen. Adding the rain at the beginning of the daily 
heating cycle makes a slight difference in the temperature, while starting the rain event 
right after the maximum heating occurred causes the overall change in depth to increase 
with time. More significant changes due to timing may have been seen if the diurnal 
heating were reduced because of significant cloud cover during the precipitation. 

mF Heavy vs. Light Precipitation (Cases 6 vs 10) 

a. Fifteen Hour Period of Precipitation at Local Noon (Case 10) 

In case 10 (Figure 12), the mixed layer temperature shows the same diurnal 


cycle with a similar 30 day change in temperature as the four hour heavy rain cycle (case 
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Figure 10. Case 8 - NPS model run using daily 4 hour precipitation of 4.5 cm/day 
starting at 0600L and | dyne/cm* wind stress. (Same labeling convention as Figure 3). 
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Figure 11. Case 9 - NPS model run using daily 4 hour precipitation of 4.5 cm/day 
starting at 1800L and 1 dyne/cm? wind stress. (Same labeling convention as Figure 3). 
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Figure 12. Case 10 - NPS model run using daily 15 hour precipitation of 4.5 cm/day 
starting at 1200L and 1 dyne/cm? wind stress. (Same labeling convention as Figure 3). 


6). The salinity decrease 1s not as sharp for each rain event because the same total 
amount of rain falls over a longer period. Over the 30 days, the salinity shows a 
decreasing trend but not as strong as the heavy rain (2 psu vice 3 psu). The lighter rain 
results in a deeper mixed layer. This accounts for the smaller salinity change for the 30 
day period because there is more mixing from below the layer. With the lighter 
precipitation, the mixed layer shoals to 9 m at the start of each new rain event but does not 
remain there throughout the rain event like the heavy rain. After the precipitation starts, 
the mixed layer deepens again. And there is an overall deepening during the 30 day period 
of about 2 m. 

Heavy vs. light precipitation does make a difference in predicting the mixed layer. 
The lighter, longer rain period slows down the total salinity decrease by allowing for more 
entrainment of saltier water from below the layer. This allows for greater deepening of the 
mixed layer over time. 

4. Wind Effects 

In these simulations, comparisons between a T value of 1.0 dynes/cm? (cases 6 and 

10) and a Tt value 0.5 dynes/cm? (case 11 and 12) are made. Reducing t by half effectively 
reduces the rate of deepening by a third. In Equation 10, the wind effects can be seen in 
the shear production terms, wind stress and entrainment shear production. 

a, Light Precipitation (Case 10 vs Case 11) 
The mixed layer temperature of case 11 (Figure 13) has a weak diurnal 

pattern like case 10 (Figure 12) but cools with time at a greater rate (1°C vs 0.7°C in 30 


days). The mixed layer salinity of case 11 follows the same trend as temperature with a 
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Figure 13. Case 11 - NPS model run using daily 15 hour precipitation of 4.5 cm/day 
starting at 1200L and 0.5 dyne/cm? wind stress. (Same labeling convention as Figure 3). 
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similar pattern but a much faster decrease over time when compared to case 10 (2.1 psu 
vs. 5.0 psu in 30 days). With t = 0.5 dynes/cm? the mixed layer depth deepens at 
approximately the same rate but shoals to 4 m instead of 9 m and has a diurnal oscillation 
of 3 to 4 minstead of 10 to 11 m. This accounts for the large change in temperature and 
salinity as this layer becomes fresher and cooler. 

b. Heavy Precipitation (Case 6 vs Case 12) 

The mixed layer of case 12 (Figure 14) has a smaller diurnal depth 
oscillation than the increased wind case 6 (Figure 8), 4 m vs 10 m daily. This 1s because 
less energy is available for deepening the mixed layer with reduced wind forcing. With the 
reduced energy for mixing, the entramment of waters from the layer directly below the 
mixed layer is reduced. This leads to the mixed layer salinity decreasing at a much greater 
rate, 7 psu vs 2 psu in 30 days. Also, the mixed layer temperature cools at a slightly 
Slower rate (0.9°C vs 1°C in 30 days) with the reduced wind. 

Wind must be taken into account, especially for light precipitation cases. For 
lighter, longer rain periods if there is no significant wind effect when the upper layer of the 
ocean is capped off by the precipitation then a very shallow, fresh, cool mixed layer is 
created. Miller's (1976) study found the same result when comparing different rain 
intensities at different wind speeds. 

C. COARE FORCING 

All cases in this section use real wind, precipitation, and radiation forcing retrieved 

from the WHOI IMET buoy at 1° 45'S and 156°E during the TOGA COARE IOP. The 


thirty day period used is 5 December 1992 to 4 January 1993. The NPS model is 
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Figure 14. Case 12 - NPS model run using daily 4 hour precipitation of 4.5 cm/day 
starting at 1200L and 0.5 dyne/cm? wind stress. (Same labeling convention as Figure 3). 


initialized with the BT temperature profile from 5 December 1992 as in the idealized 
experiments and with the climatological salinity profile from the Levitus data at this 
location. 

I Full Forcing (Case 13) 

For the first five days of the COARE case (Figure 15), there is ight wind and 
precipitation with a diurnal heating pattern. The result of this forcing was a shallow 
diurnal cycle of 1 to 5 m for the mixed layer depth. The mixed layer temperature 
increased approximately 2.5 °C while the mixed layer salinity decreased slightly. This is 
similar to what was observed in the light precipitation and light wind forcing cases earlier. 
On days 5 and 6, there 1s a large decrease in the mixed layer salinity, an increase in the 
mixed layer temperature, and the diurnal cycle is not readily seen in the mixed layer depth. 
This is explained by starting with a very shallow mixed layer and adding nominal wind and 
light precipitation, which causes the mixed layer to become insulated from the layer below. 

Increased winds and heavier precipitation events occur from days 7 to 12. The 
diurnal heating pattern is reduced by the clouds associated with the precipitation events, 
especially on days 9 and 10. With this forcing the mixed layer deepens 5 m, cools 3°C, 
and freshens nearly 0.5 psu. 

There appears to be a storm from days 15 to 17 with very high winds, large 
precipitation rates, and reduced diurnal forcing due to clouds. During this storm, the 
mixed layer temperature decreases slightly (~0.5°C) and the mixed layer deepens 5 m in 2 
days, but the mixed layer salinity remains relatively unchanged. During this event there is 


entrainment of water from the lower layer of the ocean so cooler, saltier water is being 
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Figure 15. Case 13 - NPS model run using full COARE forcing. (Same labeling 
convention as Figure 3). 


mixed into the layer. With the reduced heating, and the larger mass of water in the mixed 
layer, the effects of fresh water are mitigated. This is why the mixed layer salinity does 
not change. 

From days 25 to 29, there is light precipitation and a strong diurnal heating 
pattern, but wind mixing dominates. On day 29, the mixed layer deepens to its greatest 
depth, 19 m, during the 30 days, while the mixed layer temperature cools slightly and the 
mixed layer salinity increases due to entrainment of the underlying water. 

To check the accuracy of the NPS 1-D model, the daily model temperature profiles 
were compared with the IMET buoy temperature profiles. Looking at the upper 200 m, 
(Figure 16), the difference in the lower layers is +1°C. Because the model focuses on the 
mixed layer dynamics, there is little to no effect on the lower 150 m. Concentrating on the 
mixed layer environment, surface to 50 meters (Figure 17), one sees that the main 
difference is in the upper 3 m of the ocean. The model appears reasonably accurate, 
+1°C, for the majority of its domain (Figure 18). It does tend to overheat the mixed layer 
by 2 to 3°C when the mixed layer is very shallow and the wind and precipitation are light. 
This may seem like a large error when looking at temperature only. But when considering 
the overall heat content in the upper 3 m, the error is small. 

Di Partial Forcing - No Precipitation (Case 14) 

To provide an overall understanding as to how precipitation affects the mixed 
layer, case 14 was run using the in situ forcing from the IMET buoy except for the 
precipitation (Figure 19). The most obvious difference without the precipitation is the 


change in mixed layer salinity value, which was to be expected. If there is no 
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Figure 16. Comparison of IMET buoy and NPS model temperature profiles from 0 to 
200 m depth. Contours are temperature in degrees Celsius. 
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Figure 17. IMET buoy and NPS model temperature profiles from 0 to 50 m depth. 
Contours are temperature in degrees Celsius. 
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Figure 18. Difference between IMET buoy and NPS model temperature profiles from 0 
to 50 mdepth. Contours are temperature in degrees Celsius. 
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Figure 19. Case 14 - NPS model run using COARE forcing with no precipitation. (Same 
labeling convention as Figure 3). 
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addition of fresh water from an outside source, the mixed layer salinity will have small 
changes because the tropical Pacific’s salinity differences are small throughout the region. 
The other difference of note is the increased depth the mixed layer reaches without 
precipitation. Including precipitation in case 13, the maximum depth reached by the mixed 
layer in the 30 day period was 18 m. Without precipitation, the mixed layer depth is more 
than double that amount as it extends deeper than 50 m. To examine the more subtle 
changes, the case 14 profile has been subtracted from the case 13 profile (Figure 20). 
Besides the mixed layer salinity and layer depth, there is also a very small mixed layer 
temperature difference during the period. The largest temperature difference (0.5°C) 
occurs between days 4 and 8. This occurs at the time of a very shallow mixed layer with 
light wind and precipitation in the full forcing case which, as mentioned earlier, caps off 
the upper layer and allows for increased heating of the layer. 

J Averaged Forcing 

The averaged forcing was obtained by taking the mean of each of the COARE 
thirty day hourly forcing inputs: wind stress, precipitation, downward solar radiation, 
back radiation, sensible heat flux, and latent heat flux. This was done to see how the 
fluctuations of the forcing affected the 30 day output. 

a. Averaged Wind and Precipitation Forcing Only (Case 15) 
Including the fluctuations in the radiation terms so that a diurnal heating 

cycle would still be in place, the NPS model was run with averaged hourly wind stress and 


precipitation values (Figure 21). With the averaged wind stress and precipitation forcing 
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Figure 20. Difference between full COARE forcing (Case 13) and COARE forcing with 
no precipitation (Case 14). (Same labeling convention as Figure 3). 
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Figure 21. Case 15 - NPS model run using COARE forcing with average wind stress and 
precipitation values. (Same labeling convention as Figure 3). 


the mixed layer stayed between 1 mand 10 mdepth. The mixed layer temperature 
finished reasonably close to the real forcing case after 30 days, with both reaching 
approximately 29°C, but case 15 only had a maximum temperature of 32°C while case 13 
had a maximum mixed layer temperature of 35°C. The mixed layer salinity decreased 
from 34.5 psu to 33 psu over the 30 day period. Thus, the instantaneous forcing is 0.4 
psu saltier than the averaged precipitation and wind stress forcing 

b. Averaged Radiation, Wind, and Precipitation Forcing (Case 16) 

Using only the averaged COARE forcing (Figure 22), the same trends 
appear as with the partial averaged forcing case. The mixed layer temperature after 30 
days is 29.8°C, which is less like the real forcing mixed layer temperature than the partial 
forcing case. With case 16, the mixed layer salinity is 32.7 psu while the real COARE 
forcing of case 13 gave a mixed layer salinity of 33.3 psu at the end of the period. 

The primary finding from this series of experiments is that the variation of the 

forcing is important. To use the average heat, wind, and/or precipitation for a forcing 
period will not allow modeling of the mixed layer with reasonable accuracy. The temporal 


variation of the forcing is important. 
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Figure 22. Case 16 - NPS model run using averaged COARE forcing. (Same labeling 
convention as Figure 3). 
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IV. A BRIEF LOOK AT LARGE-EDDY SIMULATION, DISCUSSION AND 
RECOMMENDATIONS 

Several conclusions can be drawn from this study of the effects of salinity in 
equatorial mixed layers. Results from the NPS one-dimensional mixed layer model, forced 
with both idealized and in situ data, demonstrate that salinity can play a significant role in 
potentially changing the surface heat flux, with its effect on the mixed layer depth and 
mixed layer temperature. Precipitation stabilizes the mixed layer by creating a barrier 
layer, which slows down entrainment. Another important conclusion is that the 
fluctuations of the surface forcing play a significant role in the formation of the mixed 
layer. This is most apparent in the comparison between the average COARE forcing and 
the full COARE forcing. This also explains the difference between the four hour and 
fifteen hour rain events. Rain rate 1s an important consideration especially when varying 
wind mixing 1s involved. 

The NPS 1-D mixed layer model did a good job modeling the western equatorial 
mixed layer in only one dimension, but the mixed layer is governed by turbulence which 1s 
fundamentally a three dimensional phenomenon. Currently, there 1s a study in progress 
which is uses a Large Eddy Simulation (LES) model to evaluate the effects of 
precipitation on the three dimensional turbulence in the western equatorial Pacific Ocean. 
In Figure 23, the turbulence cells, known as Langmurr cells, are visible as pressure 
fluctuations from a LES model experiment at 2°N latitude, which has the same Coriolis 
value as the location used for the NPS 1-D mixed layer model in this paper. The depth of 


the vertical slice 1s 200 m to coincide with the 1-D model’s domain and the turbulence can 
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Figure 23. LES simulation at 2 


grid resolution was 5 m in each direction. 
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be seen reaching nearly to 100 m. Analysis of the LES model output enables statistical 
characterization of the three dimensional turbulence and the buoyancy flux. This study, 
and use of an embedded mixed-layer ocean general circulation model, modified to include 
salinity, will provide an even better overall understanding of salinity/precipitation effects in 
this region. 

Finally, since the western equatorial Pacific Ocean is referred to as the “warm 
poo!” and is believed to play a significant role in the global climate, future studies in this 
area should include the use of a coupled atmosphere-ocean model to study the effects of 
salinity on the feedback cycles. This should provide better understanding of the western 


equatorial Pacific’s warm pool and its role in the earth’s climate. 
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